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Abstract 

This paper studies the linear stochastic partial differential equation of fractional 
orders both in time and space variables u{t,x) = Xu{t, x)W{t, x), 

where IT is a general Gaussian noise and j3 G (1/2,2), a G (0,2], The existence and 
uniqueness of the solution, the moment bounds of the solution are obtained by using the 
fundamental solutions of the corresponding deterministic counterpart represented by 
the Fox H-functions. Along the way, we obtain some new properties of the fundamental 
solutions. 


1 Introduction 


In this article we consider the following linear stochastic partial differential equation of frac¬ 
tional orders both in time and space variables: 


j u{t, x) = \u{t, x)W{t, x), 

= Uk{x), 


dk 


t > 0, X e M'^, 
0<k<\/3]-l, a; G 


( 1 . 1 ) 


t=o 
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with [3 G (1/2,2) and a G (0,2], where [/3] is the smallest integer greater than or eqnal to 
/3. Here and thronghont the paper we denote k G N. We limit our consideration 

to the above parameter ranges of f3 and a since we plan to use some particular properties 
of the corresponding Fox H-functions which will be proved only for these parameter ranges. 
Now let us give more detailed explanation on the terms appearing in the above equation. 

dP 

The fractional derivative in time = 7 ^ is understood in the Caputo sense: 


d^fit) ; = 


r(m - /3) 

d*^ 


dP 


:fit) 



{t - 


dr 


ifm — l</d<m, 
if /3 = m , 


& 


where f > 0. A = 7 7 -^ is the Laplacian with respect to spatial variables and (—A)“A 

^ oxj 


2=1 


is the fractional Laplacian. W is a zero mean Gaussian noise with the following covariance 
structure 

¥.{W{t,x)W{s,y)) = 7 (t-s)A(x-|/), 


where both (possibly generalized) functions 7 and A are assumed to be nonnegative and 
nonnegative dehnite. We denote by /i the Fourier transformation measure of A(x). Namely, 




This Fourier transform is understood in distributional sense (see Section 2). When 7 (f) = 
6o(t) and A(a;) = 5o(a^), this noise W reduces to the space-time white noise, z/ > 0 and A are 
some real valued parameters. The given initial conditions Uk{x) are assumed to be continuous 
and bounded functions. The product u{t, x)W{t, x) in equation (1.1) is the Wick one (see e.g. 
[20]). So, the equation will be understood in the Skorohod sense. Let us point out that some 
of our results can also be extended to nonlinear equation (namely, replace u{t,x)W{t,x) in 
(1.1) by a{u{t,x))W{t,x) for a Lipschitz nonlinear function a). However, we limit ourselves 
to this linear case for two reasons: One is to simplify the presentation and to better explain 
the ideas and the other one is that we want to use the chaos expansion method. 

The deterministic counterparts of equation (1.1) have received a lot of attention and are 
called anomalous diffusions. They appeared in biological physics and other helds. Equation 
(1.1) is an anomalous diffusion in a Gaussian noisy environment. More detailed motivations 
for the study of this type of equations are given in [4, 9, 17, 27]. Here, we briefly mention 
some recent results. 

When [3 G {1,2} is an integer, a = 2, the equation has been studied by many authors, 
see for example, [1, 2, 7, 5, 18]. The work by Ghen and Dalang [ 6 ] deals with the case where 
/9 = 1, a G (1,2]. Khoshnevisan and Foondun [13] and Song [34] study a similar equation 
with the a-stable generator (—A)"A replaced by a general Levy generator. 

When /3 G (0,1), a = 2, A is replaced by a general elliptic operator, and IF is a fractional 
noise, the equation was studied in [17]. 
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When /9 G (0,1), a = 2 and W is a fractional noise, the smoothed eqnation 




u{t, x)W{t, x) 


(see (3.19) for a generalization) was stndied in [4], In a series of papers [14, 26, 27], Nane 
and his coanthors stndied the case a G (0, 2], 

The case (3 G (0,1) corresponds to the slow diffnsion (subdiffnsion). For the fast diffnsion 
case (super diffusion), i.e., (3 G (1,2), there have been only a few works. The first author of 
this paper studies in [4] the smoothed equation with a = 2, d = 1 and with space-time white 
noise. The corresponding non-smoothed equation is studied recently in [8]. Both papers [4, 8] 
deal with the nonlinear equation, i.e., p{u)W with p being a Lipschitz function. 

To study equation (1.1) the important tools that we shall use are the fundamental solu¬ 
tions corresponding to its deterministic counterpart. Let us briefly describe them. There are 
two fundamental solutions 


Z{t,x) := and Y(t,x) := Ya,^^d(t,x) 

such that the solution u{t, x) to the following deterministic equation (the deterministic coun¬ 
terpart of (1.1)) 


(d^ + ^(-A)“/2) u{t,x) = f{t,x), 

gk 


dt^ 

is represented by 


M(f, x) 


= Uk{x), 


t=0 


t > 0, X G 

0 < A; < [/d] — 1, X G M'^, 


( 1 . 2 ) 


u{t,x) = Jo{t,x) + [ ds [ dy f{s,y)Y{t - s,x - y), (1.3) 

Jo JRd. 

where / is a continuous and bounded function on M_|_ x and Uk is a continuous and 
bounded function on In equation (1.3) and throughout the paper, we denote 

r/31-i 

Jo(t,x) := 

k=0 

Equations (1.3) and (1.4) say that Z and Y are the fundamental solutions corresponding 
respectively to the initial conditions and the inhomogeneous term of equation (1.2). For some 
parameter ranges of a and (3, the fundamental solutions have been studied in [12, 24, 31, 33]. 
In Section 3.1 we shall study them for all (3 G (1/2, 2) and a G (0, 2]. In particular, we shall 
obtain some new properties such as the positivity of the fundamental solutions Y and Z. 

Equation (1.3) motivates us to study the mild solution to (1.1) (see e.g. Definition 2.1 
below), namely, the solution to the following stochastic integral equation: 

u{t,x) = Jo{t,x) + [ [ Y{t-s,x-y)u{s,y)W{ds,dy). (1.5) 

Jo JRd. 


U 


r/3]_i_fc(2/)5''^(t,x - y)dy. 


(1.4) 
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As in the classical case, the above equation can be studied by using the Ito-Wiener chaos 
expansion. To this end we need to understand well the two fundamental solutions Z and Y. 
In particular, we need their nonnegativity and some heat kernel like estimates. 

The nonnegativity of some Z’s is known. However, since Y is the Riemann-Liouville 
fractional derivative of Z, its nonnegativity is a challenging problem. There have been only 
few results: As proved in Lemma 25 of [31], Y2^i3^d with (3 G (1,2) is nonnegative if and only 
if d < 3. The one dimensional case is proved in [25], namely, DtZa^jS^i, and hence is 

nonnegative either if 1 < /? < a < 2, or if a G (0,1] and {3 G (0, 2). In this paper, we shall 
show the nonnegativity of Y in the following three cases: 

rQ;G(0,2], /3 g (1/2,1), d G N, 

iaG(0,2], /lG(l,a), d = 1, (1.6) 

la = 2, /Sg( 1,2), d = 2,3. 


This includes the above mentioned results as special cases. Let us also point out that for the 
smoothed SPDE, only the fundamental solution Z is needed, which is usually more regular 
than the fundamental solution Y. 

When (3 = 1 and a = 2, to show the solution of (1.1) is square integrable, it is assumed 
in [18] and [19] that the covariance of noise satisfies the following conditions: 

(i) 7 is locally integrable; 

(ii) Dalang’s condition [ < cxd is satisfied (see also [10, 13]). 

jR<i 1 + Kr 

For the existence and uniqueness of the solution to the general equation (1.1), Dalang’s 
condition will be replaced by the following condition: 


r /i(de) 

V1 + 


(1.7) 


It is obvious that if it is formally set a = 2 and /3 = 1, then (1.7) is reduced to the usual 
Dalang’s condition. 


The remaining part of the paper is organized as follows. We first specify the noise structure 
and present the definition of the solution in Section 2. The main results are Theorem 3.5 on 
the existence and uniqueness of the mild solution and Theorem 3.11 on the moment bounds 
of the solution stated in Section 3. The proof of these two theorems are based on some 
properties of the fundamental solutions represented in terms of the Fox H-functions. These 
results themselves are of particular interest and importance. We also list them as Theorem 3.1 
and Theorem 3.3 in Section 3. The properties of the fundamental solutions (Theorem 3.1) are 
proved in Section 4.1 by using the Fox H-functions. In Section 4.2, we obtain an expression 
of the density function for the d-dimensional spherically symmetric a-stable distribution - an 
auxiliary result (Theorem 3.3) which is used in the proof of Theorem 3.1. The existence and 
uniqueness result (Theorem 3.5) of the solution to (1.1) is proved in Section 5. In Section 6, 
we prove the explicit moment bounds when A is the Riesz kernel. 


4 




Our main results (Theorem 3.5) assume that the fundamental solutions are nonnegative. 
However, when 1 < (3 < 2 and when the dimension is high, the nonnegativity of the fun¬ 
damental solution Y is not known yet. In this case, we shall show in Theorem 3.14 the 
existence and uniqueness of the solution of (1.1) for some specihc Gaussian noise whose co- 
variance function A is the Riesz kernel. Finally, in Appendix A we collect some knowledge 
on the Fox H-function which we need in this paper. 

2 Preliminary 

Let us start by introducing some basic notions on Fourier transforms. The space of real¬ 
valued inhnitely differentiable functions on with compact support is denoted by 'D(]R'^) 
or V. The space of Schwartz functions is denoted by 5(M'^) or S. Its dual, the space of 
tempered distributions, is denoted by 5'(M'^) or S'. The Fourier transform is defined with 
the normalization 

3Fu{^) = f e~'^''''u{x)dx, 

so that the inverse Fourier transform is given by J^“^m(^) = {27i)~'^J^u{—Q. 

Similarly to [18], on a complete probability space (G, J^, P) we consider a Gaussian noise 
W encoded by a centered Gaussian family {IF(<^); (p G 'D(M+ x M^)}, whose covariance 
structure is given by 

E (IF((p) hF('^)) = / i.p{s,x)'rp{t,y)'^{s — t)k{x — y)dxdydsdt, (2.1) 

iR^xR2d 


where 7 : M —IR+ and A : —)• IR+ are nonnegative definite functions and the Fourier 

transform PK = y such that /i(d^) is a tempered measure, that is, there is an integer m > 1 
such that /]gd(l + |^P)~™'A(dO < oo- Throughout the paper, we assume that 7 is locally 
integrable and we denote 


Ct := 2 / 7 (s)ds, t > 0 . 


( 2 , 2 ) 


Let P be the completion of ©(M^. x M'^) endowed with the inner product 




' R^ X R2d 


ip{s, x)'ip(t, y)'y{s — t)A{x — y) dxd^/dsdf 


(27r) 


dsd«, 


(2.3) 


where p(p refers to the Fourier transform with respect to the space variable only. The 
mapping ip —)■ W{ip) defined on P(M+ x M*^) can be extended to a linear isometry between 
T-L and the Gaussian space spanned by W. We will denote this isometry by 


IF(0) = 


0 (t, x)hF(dt, dx), for 0 e "H. 
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Notice that if (f) and -0 are in Ti, then E (Ph(0)hh('^)) = (0,'0)-^. 

We will denote by D the derivative operator in the sense of Malliavin calculus. That is, 
if F is a smooth and cylindrical random variable of the form 




with (j)i ^ 71, f ^ C'“(M"') (namely / and all its partial derivatives have polynomial growth), 
then DF is the H-valued random variable dehned by 


DF 


j=l ^ J 


The operator D is closable from L7{VL) into L‘^{Q-,'H) and we dehne the Sobolev space 
as the closure of the space of smooth and cylindrical random variables under the norm 

||F||,,2=yEm + E|||£>F||y. 

We denote by 6 the adjoint of the derivative operator given by the duality formula 


E{6{u)F)=E{{DF,u)n), (2.4) 

for all F G and any element u G in the domain of 6. The operator 6 is also 

called the Skorohod integral because in the case of the Brownian motion, it coincides with an 
extension of the Ito integral introduced by Skorohod. We refer to Nualart [28] for a detailed 
account of the Malliavin calculus with respect to a Gaussian process. 

With the Skorohod integral introduced, the dehnition of the solution to equation (1.1) 
can be stated as follows. 


Definition 2.1. Let Z and Y be the fundamental solutions dehned by (1.2) and (1.3). An 
adapted random held {u = u(t,x) : t >0,x E such that E [u^(t,x)] < +oo for all (t,x) 
is a mild solution to (1.1), if for all {t,x) G M+ x the process 

{Y{t- s,x- y)u{s,y)l[o,t]{s) : s>0, y e 


is Skorohod integrable (see (2.4)), and u satishes 


u{t,x) = Jo{t,x) + 



Y{t- s,x- y)u{s, y)W{ds, d?/) 


(2.5) 


almost surely for all {t,x) G M+ x where Jo{t,x) is dehned by (1.4). 


The main ingredient in proving the existence and uniqueness of the solution is the Wiener 
chaos expansion, to which we now turn. 

Suppose that u = {u(t,x);t > 0,x G is a square integrable solution to equation 
(2.5). Then for all hxed {t,x) the random variable u{t,x) admits the following Wiener chaos 
expansion 

OO 

u{t,x) , (2.6) 

n=0 
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where for each (f,x), /„(■,•,t,a;) is a symmetric element in 'H®”'. Then, as in [15, 16, 17], 
to show the existence and uniqueness of the solution it suffices to show that for all (f, x) we 
have 

OO 

^n!||/„(-,-,f,a;)||^®„ < oo. (2.7) 

For technical reasons (see (2.8) below), we will assume, throughout the paper, the follow¬ 
ing properties on A: 

• A : —)■ [0, cxd] is a continuous function, where [0, cxd] is the usual one-point compact- 

ihcation of [0, oo). 

• A(x) < oo if and only if x 7 ^ 0 or -F(A) G L°°(M'^) and A(x) < 00 when x 7 ^ 0. 

With these two assumptions, according to Lemma 5.6 in [21], for any Borel probability 
measures ui^dx) and h' 2 {dx), the following identity holds, 

/ / A(x - ?/)i/i(dx)z/ 2 (d|/) =-^ /" • (2-8) 

jRd J^d [271)°' J^d 

In particular, the above result can be applied to the case when h'i[dx) = /i(x)dx and z/ 2 (dx) = 
/ 2 (x)dx for two nonnegative functions fi and /2 G L^(]R‘^). 


3 Main results 


3.1 Fundamental solutions: formulas and nonnegativity 

Our hrst result is concerned with the fundamental solutions to (1.2) stated in the following 
theorem. We need the two parameter Mittag-Leffler function Ea^g[z)\ 

OO Yi 

Eo.A^)-.= Y, v( 3?(«)>0,/?gC, zgC, (3.1) 

n=0 


where 5R(a) is the real part of the complex number a. When /3 = 1, we also write Ea{z) := 
Ea,i{z). The TT-functions appearing in the following theorem and their properties are given 
in the appendix. 

Theorem 3.1. The fundamental solutions to (1.2) are given by 


Z[t, x) 


Ea,l3,d(fy ^) ^ 


-d/2^\g] 


-1 


X 


-dRf^ 


2,3 


X 


2 “ ^vth 


(i,i),(r/3i,/3) 

(d/2,a/2), (1,1), (l,a/2) 


and 


y(t,x) 




= TT-^/^ 


X 


-V-'Hll 


X 




(l,l),(/9,/3) 

(d/2,a/2), (1,1), (l,a/2) 


(3.2) 


(3.3) 
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If P ^ (1, 2), then 



|a:|" (1,1), (1,^) 

(d/2,a/2), (1,1), [1,a/2) 


(3.4) 


The Fourier transforms of the fundamental solutions are given by the following: 


TZ{t, ■){i)= (-2-Vi'’|e 1“), 

TY(t, ■)(()= 

TZ’(t, .)(0 = Ef(-2-Utf\i\‘‘), !/ /3 e (1.2); 


(3.5) 

(3.6) 

(3.7) 


Moreover, we have the following results on the positivity of the fundamental solutions. 

(a) If (3 G {0,1], d E N and a G (0,2], then both Z{t,x) and Y{t,x) are nonnegative; 

(b) If (3 E (1,2), d E {2,3}, and a = 2, then both Z{t,x) and Y{t,x) are nonnegative; 

(c) If (3 E (1,2), d=l and a E [(3,2], then all Z{t,x), Y{t,x) and Z*{t,x) are nonnegative. 
The proof of this theorem is given in Section 4.1. 

Remark 3.2. Here are some known special cases: 

(1) When a = 2 and (3 E (0,1), it is proved in [24, 33] and in [12], respectively, that 



(3.8) 


and 



(3.9) 


which correspond to our Z 2 ,is,d{t,x) and Y 2 ,is,d{t,x), respectively. The equivalence is clear 
by applying Property A.3. For Z 2 ,is,d) see also [23, Chapter 6]. 

(2) When a = 2 and (3 E (0, 2), it is proved in [31] that 



(3.10) 


which corresponds to our Y 2 ,ii,d with u = 2. 

(3) In [25], the fundamental solution Z* i^ j^{t,x) has been studied for all a, (3 E (0,2) and 
d = 1. From the Mellin-Barnes integral representation (6.6) of [25], we see that the 


reduced Green function of [25] can be expressed by using the Fox H-function: 



(3.11) 
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where a and (3 have the same meaning as in this paper and 9 is the skewness: 16^1 < 
min(a;,2 — a). For the symmetric a-stable case, i.e., 6 = 0, this expression can be 
simplihed by using the dehnition of the Fox H-function and the fact that (see, e.g., [29, 
5.5.5]) 

= ^2“r(l/2 + as/2). (3,12) 

f (1 + asjZ) y'TT 

Hence, 




(l/2,«/2),(l,l),(l,a/2)y ’ 


x e M. 


(3.13) 


This implies that the fundamental solution in [25, (1.3)] 






(1/2,a/2), (1,1), (l,a/2) 


corresponds to our x) with u = 2. 

The proof of the nonnegativity part in Theorem 3.1 requires a representation of the 
spherically symmetric a-stable distribution from the Fox H-function, which is of interest by 
itself. The one-dimensional case can be found in [25]; see Remark 3.4 below. 

Theorem 3.3. Let X be a centered, d-dimensional spherically symmetric a-stable random 
variable with a G (0, 2]. Then the characteristic function and the density of X are, respec¬ 
tively. 


= exp (-If |“), {6 R"*, 


(3.14) 


and 


Pa,d{x)=n ^Hlj({\x\/2) 


(1,1) 

(d/2,a/2), (l,a/2) 

The proof of this theorem is given in Section 4.2. 


X G 


(3.15) 


Remark 3.4. When d = 1, the formula (3.15) yields a result in [25]. In particular, as proved 
in [25] (see (3.11)), when d = 1, we have 


Pa,i{x) 


\x 


-1 rrl,l 
-" 2,2 



(l,l),(l,a/2)\ 
(l,a), (l,a/2)) 


.Tr-1/2 

= 71 ' 


X 


-1 0 - 1,1 
-" 1,2 



( 1 , 1 ) \ 
(l/2,a/2), (l,a/2)) ’ 


where the second equality is due to (3.12) and the dehnition of the Fox H-function. 


3.2 Existence and uniqueness of solutions to the SPDE 

The following is one of the main theorems of the paper. 

Theorem 3.5. Assume the following conditions: 
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(1) Ya^j 3 ^dit,x) is nonnegative; 

(2) 13 e (1/2,2) andae (0,2]; 

(3) j is locally integrable; 

(4) in satisfies Dalang’s condition (1.7); 

(5) The initial conditions are such that for all t > 0, 

Ct := sup |Jo(s,2/)| < +CXD. (3.16) 

yeRd,, se[o,t] 


Then relation (2.7) holds for each {t,x). Conseguently, eguation (1-1) admits a unigue mild 
solution in the sense of Definition 2.1. 

The proof of this theorem is given in Section 5. 

Remark 3.6. From Theorem 3.1, it follows that the three cases in (1.6) satisfy the above 
assumptions (1) and (2). Moreover, if Uk G L°°(]R'^) for the first two cases in (1.6), or if 
Mo G L“(M‘^) and ufix) = Mi is a constant for the last case in (1.6), then by Lemma 6.1 
below, 

\Jo{t,x)\ < I |mo| 

Hence the assumption (5) is also satisfied. The Dalang condition (1.7) imposes a further 
restriction on the possible values of {a, (3) due to the spatial correlation A(a:). 

Remark 3.7 (Space-time white noise case). When the noise hF is a space-time white noise, 
i.e., jit) = 6o(t) and A{x) = So{x), then Dalang’s condition (1.7) becomes 

1 + i < 2. (3.17) 

This condition implies that (3 > 1/2. In particular, if a = 2 and <7=1, then (3.17) reduces 
to 

(3 >2/3, 

which recovers the condition in [8] and [9, Section 5.2]. 11/3 = 1 and <7 = 1, then this condition 
becomes 


a > 1, 


(3.18) 


which recovers the condition in [6]. 
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3.3 The smoothed equation 

The methodology used in the proof of Theorem 3.5 can also be used to study the following 
equation 




u{t, x)W{t, x) 


(3.19) 


with the same initial conditions as (1-1). Here jf is the Riemann-Liouville fractional integral 
of order /3 (with an abuse of the notation (3): 


I? fit) = 




{t — s)^ ^f{s)ds, for f > 0 and /3 > 0 . 


Due to the fractional integral in equation (3.19) which plays a smoothing role, the mild 
formulation for the solution can be expressed by using Z{t,x) only, namely. 


u{t,x) = Jo{t,x) + 



Z{t - s,x- y)u{s, y)W{ds, dy). 


(3.20) 


Then, using the same procedure as in the proof of Theorem 3.5, we have the following result. 

Theorem 3.8. Assume the conditions (3) and (5) in Theorem 3.5 and the other conditions 
are replaced by the following: 

(!’) Za,i 3 ^dit,x) is nonnegative; 

(2’) (3 e (1/2,1] U (3/2, 2) and a G (0, 2]; 

(f’) y satisfies 




1 + 


< oo. 


(3.21) 


Then relation (2.7) holds for each {t,x). Conseguently, the smoothed eguation (3.19) admits 
a unigue mild solution in the sense of Definition 2.1 with Y replaced by Z. 


Remark 3.9. The condition (1’) is different and is usually easier to verify than the condition 
(1) in Theorem 3.5. When (3 G (1/2,1], the condition (3.21) becomes < cxd which 

is also weaker than (1.7) (since (3 < 1). When (3 G (3/2,2), the condition (3.21) becomes 
/iRd < oo which is also weaker than (1.7) (since (3 < 2). 

The proof is essentially the same as that for Theorem 3.5, the only change in the proof 
worthy to be pointed out is that instead of computing the integral 

poo 

Jo 

we now need to compute the integral 





{-2-^uw^\^\‘^)dw 


(2/z/) (21/31-1)/^ 

|^|a(2r/31-l)//3 



1 1 





{—h's)ds 
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c 

|^|a(2r/31-l)//3 • 


(3.22) 


The integrability condition of the above eqnation at zero and at inhnity implies that /3 > 0 
and (3 G (1/2,1] U (3/2,2] (which is eqnivalent to f/?] < /3 + 1/2), respectively. Note that 
this condition on (3 is more restrictive than the condition (3 G (0,2) in [4]. 

Remark 3.10 (Space-time white noise case). When the noise hh is a space-time white 
noise(namely = d^), then Dalang’s condition (3.21) becomes 

d<a{2\(3]-l)/f3 or - + 1 < (3.23) 

a [3 p 

In particnlar, if a = 2 and d = 1, then this condition reduces io (3 <2. (3 = 1 and d = 1, 

then this condition becomes (3.18), which recovers the condition in [6]. 

3.4 Moment bounds 

In this subsection we give some upper bounds for the p-th moment and the lower bound of 
the second moment of the solution for some specihc choice of the covariance kernel. 

Theorem 3.11. Assume the following conditions: 

(1) The initial conditions satisfy condition (5) of Theorem 3.5; 

(2) {a,(3,d) satisfies one of the three conditions in (1.6); 

(3) A{x) = a; G with 


0 < K < min(2Q; — a/(3, d). 

Then the solution u{t,x) to (1.1) satisfies that for all p > 1, 

E [\u{t,x)f] < C^dfexp (f(aC'tCC'*(2/z/)'^/“(27r)-^)2<«/3-“-/^-p2'"^^ , (3.24) 

where Ct and Ct are defined in (2.2) and (3.16), respectively, C = C{a,/3,K) > 0, and 

a = Ti2(3-l-(3K/a) and d=[ El^i-\^n\^r^d^, 

and Ck appears in the Fourier transform of \x\~'^, i.e., p{df) = 

In particular, if ■j is the Dirac delta function and if the initial data uo{x) = uq > 0 is a 
constant and ui = 0 when (3 > 1, then for some constant c = c(a, f3, k) > 0, 

E [\u{t,x)\^] > eng exp (^t (2/z/)'^/“)2/3-iV/-j . (3.25) 

The proof of this theorem is given in Section 6. The same method can be used to obtain 
the moment bound for the solution to the smoothed equation (3.19). 
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Remark 3.12. When (3 = 1 and a = 2, equation (1.1) is reduced to the multiplicative 
stochastic heat equation (1.1) considered in [18]. In this case the exponent of p in (3.24) 
becomes 

2af3 — (3 k 4 — k 

2a(3 — a — 13k 2 — k ’ 

which is the same as in [18, Theorem 6.1, inequality (6.1)] (with k = a). If we assume 
7 (f) = t~^, then Ct = The exponent of t in (3.24) is 


1 + (—/3 + 1 ) 


a 


2a(3 — a — (3 k 


4:-2(3-k 
2- K ’ 


which is the same exponent of t as in [18], inequality (6.1). Hence, we conjecture that the 
bound (3.24) is sharp. 

Theorem 3.13. Assume the conditions (1) and (2) of Theorem 3.11, and assume 
(3’) A{x) = X elSA with 

0 < K < mm{a/(3, d). 


Then the solution u(t,x) to the smoothed equation (3.19) satisfies that for all p > 1, 

E [u{t, xY] < C^Cfexp it [aCtC'C'#(2/z/)'^/"] , (3.26) 

where C = C{a,/3, k) > 0, Ct is defined in (3.16), 


C# = r(2 r/?l - 1 - /?K/a) and C=[ Bftr/il (" 


and is as defined in Theorem 3.11. In particular, if 'y is the Dirac delta function and if 
the initial data Uo{x) = Uq > 0 is a constant and Ui = 0 when (3 > 1, then for some constant 
c = c{a, (3, k) > 0, 

E[\u{t,x)\^] >culexp(^t ■ (3-27) 

The proof of this theorem is a line-by-line change of the proof of Theorem 3.11, and we 
leave it to the interested reader. 


3.5 Case 1 < (3 <2 and d>2 

When 1 < (3 < 2 and a 7 ^ 2, we could not show the nonnegativity of Y{t,x) for high 
dimension [d > 2) (see Theorem 3.1 (b)). However, with a slightly different approach, it is 
possible to obtain similar results to Theorem 3.5 for the Riesz kernel case. 

Theorem 3.14. Assume the conditions (2), (3) and (5) of Theorem 3.5, and assume 

(4’) A{x) = \x\~'^, X with 


0 < K < min(2a — a/(3, d). 
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Then relation (2.7) holds for each {t,x). Consequently, equation (1-1) admits a unique mild 
solution in the sense of Definition 2.1. 

This theorem is proved in Section 7. 

Remark 3.15. It is easy to see that the condition A(a;) = \x\~'^ with D < k < 2a — a/fi 
implies Dalang’s condition (1.7). Condition k < disto guarantee that A is a positive function. 


4 Fox H-functions: Some proofs 

4.1 Proof of Theorem 3.1 

The proof of Theorem 3.1 will be based on following lemmas. 

Lemma 4.1. The function Za^i 3 ^d{t,x) has the Fourier transform given by (3.5). 

Proof. The proof needs relation (A. 19) between the Mittag-Leffler function and the Fox H- 
function. We hrst note some special cases. The case where /3 G (0,1], a = 2 and d E N can 
be found in [24, Section 4] or [33]. For fi G (0,1] and for general a, one can simply replace 
I'CP by I'C|“ ill the argument of [24, Section 4] and then use (A.19) to obtain (3.5). The case 
where d = 1, /? G (0, 2), and a G (0, 2) is proved by [25]. For the general case, solves 

+ ^(—A)“/^j u{t,x) = 0, f > 0, X G M, 

m( 0, x) = (5o(x), ifmG(0,1), 

m( 0 , x) = 0 , f M(t, x) = ( 5 o(x), if m G [ 1 , 2 ). 

Hence, the Fourier transform of Z^^y^d satishes 


di’Tzit.-m 




dm 

dt'^ 


-Fz{t,-m 


i=0 


1 . 


This equation can be solved explicitly (see, e.g., [11, Theorem 7.2, on p. 135]) as 


FZ{t,■)((,) ^ ITEy-vmyi), 


which gives immediately (3.5) when m = 0. When m = 1, the integral can be evaluated by 
[30, (1.99)] to give 

DZ{fi-m = tEy^fi-v\i\Hf^/2). 

This completes the proof of Lemma 4.1. □ 

Lemma 4.2. The function Za^y^d{t,x) can be expressed in (3.2). 

Proof. Following Lemma 4.1, we need to compute the inverse Fourier transform of (3.5). 
Instead of hnding the inverse Fourier transform, it turns out that it is easier to verify that 
the Fourier transform of (3.2) is equal to the right hand side of (3.5). Let now Z be dehned 
by (3.2). 
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Case I (i = 1. Notice that x h-)■ Z{t,x) is an even fnnction. We have that 


TZ{t, ■){0 = 27r-^/2^r/31-i / ^ 3 . 


X 


^ I 2^-^vtP 


(i,i),(r/3h/3) 

(1/2,a/2), (1,1), (1,a/2) 


cos(a;^). 


Write the cos(-) fnnction in the Fox H-fnnction form by (A.20) and then apply Property A.5 
to both Fox H-fnnctions and Property A.4 to the Fox H-fnnction coming from cos(-): 


.In X \ 


X 


2l-l/ajyl/Q^/3/a 


(1.1/A. (r/31,/3/a) 

(l/2,l/2), (1,1/a), (1,1/2) 


X h: 


0,1 

2,0 




(1,1/2), (l/2,l/2) 


Now we will apply Theorem A. 8 . Notice that both condition (1): = {2 — (3)/a > 0, = 0, 

A 2 = —1 7 ^ 0 , and condition (A.16): 

Ai = 0 , A 2 = 00 , Bi = min(l, a), B 2 = 0 , 

of Theorem A .8 hold. Hence, Theorem A .8 implies that 




(1,1/a), (1,1/2), (l/2,l/2),(r/31 ,/3/a) 
(l/2,l/2), (1,1/a), (1,1/2) 


, -1 


2-l/ajyl/a^/3/a 

Then apply Properties A.5 and A.3 to simplify the above expression: 

(H/fr) 

This proves the lemma when d = 1. 


(l,l),(l,a/2),(l/2,a/2),(r/31,/3) 
(l/2,a/2),(l,l),(l,a/2) 

(0.1) 

(0,l),(l-r/31,/3) 


Case II d > 2. Becanse the fnnction x h->• Za^i 3 ^d{t,x) is a radial fnnction, it is known that 
(see, e.g., [36, Theorem 3.3 on p. 155]), 


POO 

TZ{t,■)((,) = / AxHlj 

Jo 


/ 


(1.1). (1/31 ./3) 

(l/2,a/2),(l,l),(l,a/2) 




where Jv{x) is the Bessel fnnction of the first kind. Then we can apply Theorem A. 10. 
Similar to the previons case, all conditions are satisfied with the condition a* = 2 — /3 > 0. 
Hence, 


(ixtfKr) 


-1 (l,l),(l,a/2),(d/2,a/2),(r/31,/3) 

(d/2,a/2), (1,1), (l,a/2) 

(0.1) 

( 0 , 1 ), ( 1 -[A,/3) 


where the second eqnality is dne to Property A.3. This completes the proof of Lemma 4.2. □ 


Note that the case d>2 can be also proved by writing the Bessel fnnction in the Fox H- 
fnnction form by (A. 21 ) and then applying Theorem A. 8 . We leave the details for interested 
readers. 
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Lemma 4.3. The fundamental solutions Ya^p^it^x) and Z^ i^ j^{t,x) are given by (3.3) and 
(3.4), respectively. 


Proof. We first prove the expression for By Section 2 of [12], we know that Ya^y^dif^x) 

is the Riemann-Liouville fractional derivative in t of Za,j 34 {t,x) of order \j3'\ — /3. Notice that 
Za,p,d{.^,x) = 0 for |a:| 7 ^ 0. Denote the Riemann-Liouville derivative of order fd G (0,1) by 
i.e., 

(</) W = Sf(T^/ 

By Property A.4, 


Zo,^y^dit,x) = Tl ^\x\ 


-)a—l 


Ut^ 


\x\ 


(l-d/2,a/2), (0,1), (0 ,q:/ 2) 
(0,l),(l-r/31,/3) 


Because a* = (2 — /?) + (2 — a)/2 > 0, we can apply Theorem A.11, 


r)\d\—dy (± \ _ (i/2 I I —cii/3—1 rrliS 

Z/oi,l3,d\ty Xj 71 |X| t -^4,3 


\X\ 


(l-r/31,/3), (l-d/2,a/2), (0,1), (0,C(/2) 
(0,l),(l-r/31,/3),(l-/3,/3) 


Then we use Properties A.3 and A.5 to simplify the above expression to obtain (3.3). The 
expression for Z^y^ can be proved in a similar way. □ 

Lemma 4.4. The Fourier transforms of Ya^y^it^x) and Z^ y^^(t,x) are given by (3.6) and 
(3.7), respectively. 


Proof. We hrst consider Ya^y^d- From Lemma 4.1 and the proof of Lemma 4.3 it follows 


Because a* 




( 0 , 1 ) 

(0,l),(l-r/31,/3) 


2 — (3 > 0, we can apply Theorem A. 11 to obtain that 


■)(«) = (2-vt»iei 


(l-r(31,/3),(0,l) 

(0,l),(l-r/31,/3),(l-/i,/3) 


This is simplihed to (3.6) by Properties A.3 and A.5. The identity (3.7) can be obtained in 
a similar way. □ 


Lemma 4.5. For all p. > 0 and 0 < 6 < min(l,/i), the following H-function is nonnegative: 

^ 1 ))-°’ for all xe^. (4.2) 

Proof. We only need to prove that the following function is nonnegative 


f{x) = |a;| ^hI’^ (^|a;| 


( 1 . 1 ) 


X G 


By Theorem A.9 and equation (A. 19), the Laplace transform of / is equal to 


dxe ^""/(x) = Ee,f,{-z). 

By [32], we know that the above Mittag-Leffler function Ea^y{—z) is completely monotonic 
if and only if 0 < a < min(/?, 1). Then the Bernstein theorem (see, e.g., [38, Theorem 12a]) 
implies that the function f{x) is nonnegative. □ 
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Lemma 4.6. The nonnegative statements in Theorem 3.1 hold true. 

Proof. We first prove the case (a). In this case, fd G (0,1]. When /3 = 1, Z and Y coincide and 
they are alpha stable densities. Hence, we need only consider the case fd G (0,1). Because 
linij^o 3;) = 0 for all |x| ^ 0, and noticing that = Id (see, e.g.. Theorem 

3.8 of [11]), we see that 

Za,fi,d{t,x) = ll~^d^~^Zc,,p^d{t,x) = ll~'^Ya,g,d{t,x). 


Hence, it suffices to show the nonnegativity of Ya^g^dif^x). Notice that by Property A. 5, 


J : 


ds 


1.1 

2 


X 


)a-l, 


( 1 , 1 ) 

(d/2,«/2), (l,a/2) 



-h 


(/ 3 ,/ 3 )\ 

( 1 , 1)7 



\X 


ajg 


2(a-l)//3^1//3' 


(l,l//3) \ ^1,0 (1 

(d/2,a/2/3), (l,a/2/3)J 14 


(/3,1) \ 
(l,l//3 )J ■ 


Now we check conditions in Theorem A. 8. Condition (1) of Theorem A.8 is satished because 
al = l/(d > 0 and = {I — (d)/(d > 0. Condition (A. 16) of Theorem A.8 holds because 


Ai = 0, A 2 = (d, Bi = djd/a, B 2 = 00 . 


Hence, Theorem A.8 implies that 


J = (d-^H. 


2,1 

2,3 


X 


a//3 




(l,l//3), (/3,1) 

(d/2,a/2/3),(l,l//3),(l,a/2/3) 


Then by Property A. 5, we see that 


poo 

Y{t,x) = /dTr~‘^^‘^t^~^\x\~‘^ / dss~^Hl’l 

Jo 


X 


> Q :—1 


U 


( 1 , 1 ) 

(d/ 2 , 0 / 2 ), (l,a/2) 



(/3,/3)\ 
( 1 , 1)7 ■ 


By Lemma 4.5, the second H-function in the above equation is nonnegative. On the other 
hand. Theorem 3.3 tells us that the hrst H-function is nonnegative. Thus, Y{t,x) is nonneg¬ 
ative. 

As for the case (b), it is known from [31] that Y 2 ^i 3 ^d is nonnegative for d < 3. By the 
same argument as in the proof of (a), .^ 2 ,; 9 ,d is also nonnegative. 

Finally, for the case (c), it is proved in [25] that x) is nonnegative. By the same 

reason as in the proof of (a), and .^a,/ 3 ,i are fractional integrals of 2'*^^ of orders I — ft 

and 1, respectively. Therefore, both Fa_^ 1 and .^a,/ 3 ,i are nonnegative as well. The proof of 
Lemma 4.6 is now complete. □ 

Proof of Theorem 3.1. The Theorem 3.1 follows from the above lemmas. □ 


4.2 Proof of Theorem 3.3 

Proof of Theorem 3.3. The characteristic function (3.14) of X is proved in [37, (7.5.3) on p. 
211]. For the density pa 4 , we need to compute the inverse Fourier transform. From [37, 
(7.5.5)] this inverse transform is 

poo 

Pa,d(r) = / e-*“ J(,_ 2 )/ 2 (rf)f''/Mf . 

Jo 
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By (A.21) and (A.22), we have that 


t(‘'+2)/V(d-2)/2(rt) = (2/r)("+2)/2^i.| 


0 I r I 


(d/2,1), (1,1) 


and 


Hence, 






(0,1/a) 


_ ^-d/2^-d I °° t-ljjhO ^ ^ 


Pa,d(r) = 71 / r 


0,2 


(d/2,a/2), (l,o/2) 


(« 


(0,1/a) 


dt. 


Application of Theorem A.8 to evalnate the above integral yields the theorem. 


□ 


5 Proof of Theorem 3.5 


Proof of Theorem 3.5. Recall that Jo(t, x) dehned by (1.4) is the solntion to the homogeneous 
equation. Using an iteration procedure as in [18], we have 

/n('Sl) 3^1) ■ ■ ■ ) Xn., t, x) Xi, ' ' ' , t, x)To(S(t(1)) ^(t(1)) 

where 

Qn^^li ^li ' ' ' 1 ^ni ^ni U 77 !^^^ ^cr(n)y ^ ^a{n)') ' ' ' '5cr(l)) ^a{2) ^(t( 1 )) ) 

and cr denotes a permutation of {1,2, ■ ■ ■ , n} such that 0 < So-{i) < ■ ■ ■ < So-(n) < t. Fix f > 0 
and X e M'^, set fn{s, y, t, x) = /n(si, i/i, • • • , s„, i/„, t, x). Then we have that 




= n\ 



n n 

dydz fnis,y,t,x)fnir,z,t,x)YlHyi “ 

i=l i=l 


where dy = dyi---dym the differentials dz, ds and dr are dehned similarly. Set fi{df) := 
ut.pm- Using the Fourier transform and Cauchy-Schwartz inequality together with 
(2.8), we obtain that 


n!||/„(-,-,t,x)||^®„ / / J^ 9 n{s,-,t,x){f)J^gnir,-,t,x){f)y{df)Y\-f{si 

(271) J[0,t]2" jR-^d 


i=l 


c? 


n\ 


- 


X 


'[0,1]2" \jR-n.d 


(J^^„(s, •,t,x)(0)V(d0 

X 1/2 


1/2 


(*A^n(r, ■,f,x)(0) h(dO) Jj7(s7-I'Odsdr, 

2 i=l 


rj)dsdr 

(5.2) 
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where the constant Ct is dehned in (3.16). Thns, thanks to the basic ineqnality ab <2 + 

6 ^) and the fact that 7 is locally integrable, we obtain 

^^2 f f ^ 

n!||/„(-,-,t,a;)||^®„ < —■,t,a;)(OlV(dO 

(^TT) J[0,4]2n jRnd 

^2 (^n I f r 

J[0,t]" JR^d, 

where the constant Ct is dehned in (2.2). Fnrthermore, from the Fonrier transform of Y{t, ■) 
we can check that 


\J^gn{r,-,t,x){0\‘^ 


/ j\2 J_ J_ ('5cr(j+l) 2 S(T(i) )'^ l'Co-(i) + ' ' ' 

^ d j=i - 

where we have set Sa-{n+i) = t. As a conseqnence, 

[ |*^^n(s, ■,^,a;)(OlV(dO 


< 


< 


< 


( n !)2 11 

^ > i=\ 

1 


JJsnp 


(n!)2 n 

^ ' i=l 


Rd 


n « 

( [(^-h+1) “ - 2"V(s<,7+i) - s,(i))^|ea(i)r)] V(dCa(*)), (5.3) 

[n\y J^d 

where we have nsed the fact that = 1 and that Y and A are nonnegative to get rid 

of the snpremnm in 7 . Therefore, nsing Fonrier transform again we have 

^2^71 r r 

n!||/„(-,-,t,x)||^®„ </ /x(dO / ds 

jRnd JT„(q , 

n ( 5 - 4 ) 

X Yl (si+i - - 2"V(si+i - Si)^|^i|") , 

-, — 1 


where T„(f) denotes the simplex 

Tn{t) := {s = (si, ■ • • , s„) : 0 < Si < ■ • • < s„ < f}. (5.5) 


By the change of variables Sj+i — Sj = tCj for 1 < i < n — 1 and t — Sn = Wn, 'vj'e see that 




(27]-)nd 


i=l 


19 
















where 


St,n = {(wi, • • • , Wn) e [0, OoY : Wi-\ -h Wn < t}. 

We take TV > 1 which will be chosen later, and let 

““d e R" : 1^1 < N}. 

J\^\>N Isl ' 

Let / be a subset of {1, 2, • • • , n} and = {1, 2, • • ■ , n}\I. Then we have 


and Dn = G : |^| < TV}. 




St,n A — 1 


Vwf|6r)(l{|c,|<7V} + l{|€.|>iv})dwi/i(dCi) 


St.n 7 = 1 


dw / /i(dOn^A/3(-2 ^l{h,|<iv} 


/C{1,2,--,n} is/ 


where dtc = dtci ■ ■ ■ dwn- For the indices i in the set I, for some constant Cjs > 1 (one may 
choose C /3 = r(/T)“^) 


4„(-2-Vufi&i») <c^. 


Now using the inclusion St,n C Sj x S'/"" with 


S'/ = G /) : tc/ > 0, '^^Wi < t| and S'/"" = G C) : tc/ > 0, < t| , 


we obtain that 


St,n i — 


Jituf ^Elf{-2 ‘!/ttjf|/i|")dTOj^(dfi) 


<C E [ f-ido 


/C{1,2,--- ,n} 


's/xSf 


Xfl^f \\^,\<N}Y[wf ' 1 {|^^.|>/ v }^^,/ 3(-2 VwflOr) 


Furthermore, one can bound the integral over S'/"" in the following way 


C./C -* 

s is/- 


n^f ior)dw^i < ^eIp{- 2 ^uw^\^j\^)dwj. 


|/C| J.J. 

+ i6/‘= 


Then make the change of variables w^\^j\°‘ —)■ Uj to obtain 
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1 


< 


csn 


j£l<^ 


|^.|2a—o//3 ’ 


where 


1 




Note that the integrability of the above quantity at zero and at inhnity implies that (3 > 1/2 
and (3 > 0, respectively. Thus we have the following bound. 


J St,n 
< 






E C';' I ]lw?-^-dw,-{,.{i€W‘:K\<N}f cl';' 

/C{l,2,-,n} i&I “'l5jl>NVie/'= Jgjc Isjl 


< 


E 


< 


r(( 2 /?-i)|/| + i) ^ ^ 

^ /„\ +(2/3-l)fc 

k=0 ^ ' 


/C{l,2,-,«} 


T{{2^-i)k + iy 


where = max(C'/ 3 ,and C]\f and are dehned in (5.6). Observing the trivial 
inequality (//) < 2**, we have 


OO OO 


n = 0 ^ U—n rr^—U \ / 


/c=0 n=k 


OO OO 


^(2/3-l)fc 


T{{2y-l)k + l) 


r^k ^n—k 

X>n'^N 


n2 _^ ^^ 

<^EE 


t(2/3-l)fc 


Choosing N sufficiently large so that 2C^CtCj^ < 1 yields 


Y\k r^—k(c\^ ^ \n 


°° Pi2 °° j.(2/3—1)A: (OC C C 

gn!||/„(-,-,f,a;)||^«„ < E y{{2(3 - l)k + _ Ic[c,Cn ^ ' 

This proves (2.7), and thus the existence and uniqueness of the solution. 


□ 


6 Proof of Theorem 3.11 

Lemma 6.1. Suppose that the initial conditions Uk{x) = Uk are constant. Then under the 
three cases of (1.6), we have that 


Jo{t,x) = 


uo if (3 E (0,1], 

uo + t^~^ui if (3 E (1,2). 


( 6 . 1 ) 
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Proof. By Theorem 3.1, we know that under the hrst two cases of (1.6), the fundamental 
solutions are nonnegative and hence. 


1 / 31-1 . 1 / 31-1 

^ / d^Z{t,x-y)(ly = ^ ^ [d^ (0), 

I —n •/R'* I —n 


which is equal to the right hand side of (6.1). As for the last case in (1.6), because Z is still 
nonnegative, the contribution by uq can be computed in the same way. However, we do not 
know whether Z* is nonnegative, and thus we cannot use the Fourier transform arguments 
to compute the contribution by Ui. Instead, we compute it directly: 



(f, x)dx 


Sd-lT^ 


-d/2 





(l,/3) 

(d/2,1), (1,1) 


da:. 


where 


Sd-i 


Wm' 


Then by Theorem A.9, we have the following Laplace transform; 


( 6 . 2 ) 


9{.z) 



e-^^x-^ 


H 


2,0 

1,2 



(l,/3) 

(d/2,1), (1,1) 


dx 



(l-d/2,l),(0,l)\ 
(0,a),(0,/3) ) 


Then by Theorem A.7, 


9(0) 


h* — 

'''10 ~ 


r(<i/ 2 ) 

2 


Putting these identities together, we have that 


/ Zlp^^{t,x)dx = l. 

dRd 

This completes the proof of Lemma 6.1. □ 

Proof of Theorem 3.11. Since A(x) = we have y.{df) = for some coefficient 

Ck] see, e.g., [35]. We begin with the upper bound. By the hypercontractivity property of 
the n-th chaos, i.e. 


\\In{fn{-,-,t,x))\\LPin) <{P- l)"||4(/n(-,-,^,a3))||L2(0). 
On the other hand, from the proof of Theorem 3.5 (see (5.4)) it follows 

\\In{fn{-r,t,x))\\l2(^^^ = nl ||/n(-,-,t,a/)|£®n 

-kps- / !!(*<+■ 

[zn) Jr„(t) JM.”d 


(6,3) 


< 




(2ir: 


nd 
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~ [Y'kY^ 

where Ct is defined in (3.16), 


h^{t) 


C-.= ElY-lvnivr^'dv = ElY-nt^ dt, 

Jr‘^ Jo 

and Sd-i is defined in (6.2). According to the property of the Mittag-Leffier function at zero 
and infinity, if 0 < k < 2a, then the above constant C is finite. Then, under the condition 
that K < a{2 — 1//3) (this condition implies 0 < k < 2a), the integration over ds can be 
evaluated explicitly; see [18, Lemma 4.5]. Hence, 

1 _ f(2/3-l-^)ri. (2luY^/°‘ 


where C* := T{2f3 — 1 — (Jk/o). Denote 


0 * := 


(27r) 


-.C^C.CtC 


Thus we obtain 


0-/2^(/3-i-|^)n 

|4(/n(-,-,t,x))|U2(^)<a- 


T((2/l-l-^)n + l )2 

This bound together with the hypercontractivity implies that 


\In{fn{'} 'jt, x))\\i^P(^Q'j 'El Ct 


r((2/?-l-^)n + l)i 


Therefore, 


\u{t,x)\\LP{n) < \\In{fn{-r,t,x))\\LP{n) <CtY 


0 / t p2 


n=0 


„^„r(2te + i)i 


(6.4) 


(6.5) 


where 


e ■= (3 -1/2- 13 k/{ 2a). 

Then by the fact that r(l + 2x) > r(l + x)"^ for x > —1, 

0 / t 

n=0 

< CCtew (t(C'«G5C'42/n)''/“(27r)-‘^)2-^--^«p2-/3-“-/3«j ^ 


k(i,i)iu.,n, < a %/!l% “ (qL^V'" 


( 6 . 6 ) 


for some positive constant C = C{a, j3, k), where in the last step, we have used the asymptotic 
property of the Mittag-Leffier function (see, e.g., [30, Theorem 1.3]). 
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Now we consider the special case when 7 is the Dirac delta fnnction. By Lemma 6.1 and 
the assnmptions on the initial conditions we have 

Jo(t, X) =Uo + = Uq. 

From the proof of Theorem 3.5 (see (5.2) and (5.4)), we see that 

1 r r 

1 / , / , ^ TT . \2/3-2 


1 r r 

l|4(/n(-,-,t,2:))||i2(o) = —/ ds / d^ JJ(s<,(i+i) - S^P)) 

X (-2“^z2(s„p+i) - S^P))^|^^P) H-h ^tT(i)r) 

[ ds [ n(si+i - SiY^- 

JTnit) JR’^d 




(27r)»^'^ 

X El^ (-2-V(s,+i - S,)^|6 + • • • + 61“) I6r“" • 

Then by the change of variable 6 + ''' 6 = hi replacing by we obtain that 


^ic: 


\\In{fn{-rAxmUn) = ji^, 


I I TT('5i+l 

Jm’-d 


2/3-2 


> 


ulc: 


(27r: 


nd 


> 


'“oC; 


(2n: 


nd 


X (-2 V(si+1 - Si)^|7,|") \7]i - irii_i\'^ "^d^idsi 
X Ej^^ (-2“V(si+i - Si)^|7,|") \T]i - 7i_i|'"“'^d^idsi 
X E^^^ (-2"V(si+i - Si)^|7i|") |7i|'"“‘^d^idsi, 


where 70 = 0. Then with another change of variable (i//2)^/"(sj+i — SiY^^rji —)■ 7 ,, and by 
the same reasoning as before, we obtain that 




||-fn(/n(-, •,t,a;))||£,2(f^) 6^27]-)nd 


/ / ^ “ ^i/3(-|h*r) IhiT '^d^idsi 




« n 

(2/1/)"^/ ]^(si+i - Si)^^"^"^dsi 

■^'^nW i=l 


(27r)^"' I 2^^ 

_r(2/3-l-^)(2/z/)W«„2Q^n('4^^-nd(7n 

“ r(n(2/? - 1 - ^) + 1) ■ 

Therefore, by the asymptotic property of the Mittag-LefHer fnnction. 


~ ng aC'(47r)-'^a) r(2/3-i-i^)(2/j,)^ 
E [u{t,xf] —h 


n=0 


r(n(2/3 - 1 - /^) + 1 ) 

>CM 2 exp (^(a(5(47r)-'^C'62/i/)'^/“)2^-i-^-/“f) , 
for some positive constant c = c{a,f3, k). This completes the proof of Theorem 3.11. □ 
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7 Proof of Theorem 3.14 


In this section, C = Ca^/s^... denotes a positive constant, possibly dependent on a, /3, d, z/, • • •. 

Lemma 7.1. Assume that /3 G (0,2), a > 0 and d G N. Then there is a nonnegative eonstant 
Ca,i 3 ,d such that for all 0 < ( < niin(d/a, 2), 


-^ 2,3 


(1,1), (A/3) 

(d/2,a/2), (1,1), (l,a/2) 




for all z >0. 


Proof. We first note that condition (A. 6 ) is satished. Because a* = 2 — (3 > 0, we can apply 
Theorem A .6 to obtain that 


H. 


2,1 


2,3 


(1,1), (/3,A 

(d/2,a/2), (1,1), (l,a/2) 


=0{l/z); 


z ^ oo. 


As for small z, note that the poles of r(l + s) are 


A:={-(l + fc):fc = 0,l,2,---}; 


and those of r(d/2 + q;s/2) are 


I a 



To hnd the leading term when ^ 0, we need to hnd the hrst nonvanishing residue of 

TL 2 '^{s)z~^ at poles AU B, where 



r(d/2 + as/2)r(i + s)r(-s) 

T{(3 + ^s)T{-as/2) 


(7.1) 


Case I. When d ^ a and d ^ 2a, then the leading pole {I = 0) in i? does not coincide 
with the hrst two poles {k = 0, 1) of A. Hence, the asymptotic expansion in Theorem A.7 
( 1 ) implies that 


H. 


2,1 

2,3 


(1,1). (AA 

(d/2,a/2),(l,l),(l,a/2) 


+ hl^z + hl,z^ + 0{z^) + 0{z^/^)- 


z ^ 0. 


Here h*- are dehned in (A. 12). Notice that = 0 due to the presence of the parameter 
(/S,/3). Hence, 


z:r2,l 

-^ 2,3 


(1,1), (/3,A 

(d/2,a/2),(l,l),(l,a/2) 


= 0 (^ 2 ) + 


2 ; -)■ 0 . 


Case II. Now we consider the case when d = a. The hrst pole in A U H is s = —1, which 
is of order 2. Now we compute the mentioned residue: 


Res [^^ 3 ( 5 ) 2 ;*] = lim [{,s + lf'H^fl{s)z "]' 

5 = —1 ’ S—¥ — l 

= lim \{s + lY'Hl{s)'H 2 i^)z~^y 

S—>-l 
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= lim [Ul{s)'ni{s) + Ul{s)ni{s)' - n\{s)Ul{s) log^], 

s^—1 

where is defined in (7.1) and 

H*(s) = (s + l)2r((l + s)d/2)r(l + s) and Hl{s) = 

Now simple calculations show that 


r(-g) 
r (/3 + /3s)r(-ds/2)' 


'Hi(-l) = lim 'H\{s) = 3 = lim — --r = 3, 

s-s>-i d s-s--! ((1 + s)(i/2)(l + s) d 

= lim n* 2 {s) = 0, 

s—> —1 




= lim 


r(-s) 


s=-l 


-ir(-ds/2) Vr(/ 3 (i + s)) 

r(i) i^m + s)) 


T{d/2)s^-i T{/3{l + s)) 

/g 

nd/2) ’ 

where 'ip{z) is the digamma function and the last limit is due to (5.7.6) and (5.7.1) of [29]. 
Thus, 


ni{-l)H;{-l) = Hli-iyuii-l) = 0 and 77*(-1 )-h;(- 1)'= 


2/? 




dr(d/2) r(i + d/2)‘ 


Hence, 


/9 2; 


Res \'H 2 \{s)z^) = 7/ \ - 

s=-i L 2,3V ; j r(l + d/ 2 ) 

Therefore, by the definition of the Fox H-function, 

2 , 1 / (1,1), (/3,/3) \ n( \ 

^ - J =0[z)- z^O. 


mi I z 


(1,1), (/3,/3) 

((i/2,d/2), (1,1), (l,d/2)y 


Case III. As for the case d = 2a, the first pole in A U i? is s = —1, which is a simple 
pole. As calculated before, the residue at this pole is vanishing, /i^qZ = 0. Hence, we need to 
consider the next pole at s = —2, which is a pole of order 2. Use the asymptotic expansion 
(A. 11) to obtain that 


m'3 I 


(1,1), (/3,/3) 

((i/ 2 , 0 / 2 ), ( 1 , 1 ), (l,a/ 2 ) 


) =0(z^ logz); 


z ^ 0. 


Finally, because a* > 0, by Theorem A.2, our H-function is a continuous function for 2; > 0. 
With this, we complete the proof of Lemma 7.1. □ 
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Lemma 7.2. For all a G (0, 2], d G N and n < min{2Q;, c?}, one can find ( < min((i/Q;, 2) 
and a nonnegative constant C (independent of a) such that 


where 


\x — a\ ^Q{x)dx < C < oo for all a G 


©(x) = 


\xfi+d |a;|rf-Ca 


Proof. We divide the integral domain into {|a;| < 1} and {|a;| > 1}. Over the domain 
{|a;| < 1}, we have 


/ \x — a\ '^0(a;)da; < / \x — a\ — 


d—(a 


dx 


\x — a 


J |a:|<l, |a:|<|a:—a| 

< / lx 

1 




^- dx f It — o\ ^ - 

N<l,N<|x-a|' ■ ' |x-a|<^-C« 


dx 


< 2 


l\z\<l 


I ^ I K+d—(a 


dx < C. 


The last inequality is valid since we can choose ( sufficiently close to m.m.{d/a,2) so that 
n + d — (a < d. On the other hand, over the domain {|t| > 1}, we have 


/ |x —a| ''0(x)dx < / |x —a| —-jdx 

'n>i ~ J\.\>i 

< [ |x - *—dx + [ Ix-al""" , 

yix-a|>|x|>i yix|>ix-a|>i 


dx 


+ / |x — a 

J |a:|>l>|x—a| 

< 2 [ + [ 

~ J\Z\>1 J\Z\<1 


;dx 


xfi+d 

|x|-Mx < C. 


Note that the above constant C does not depend on a. 


□ 


Lemma 7.3. Assume k < min{2Q;, d}. Then for all s,r > 0 and X 2 , 1/2 G 1^^, we have that 
[ \Y{s,xi-X 2 )Y{r,yi-y 2 )\ \xi - yi\~'^dxidyi < Ca,g,d,u,K {s rf, 


where C does not depend on X 2 and y 2 G and 




2a 
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Proof. We use the notation Q{x) in the previous lemma. By Lemma 7.1 and the expression 
of Y through Fox H-function (3.3), we see that for any ( < mm{d/a,2), there is a constant 
Ca,i3,d,u,c, such that 


\Y{t,x)\ 


< 




I tl^/a \ 


I X laC+g I 1 
I t/3/c I W i 





Therefore, by Lemma 7.3, we have 


\Y{s,Xi-X 2 )Y{r,yi-y 2 )\ \xi-yi\ ""(Ixidyi 






jR' 

jRd 

By symmetry, we also have 

|F(s,xi -X 2 )Y{r,yi - y 2 )\ ki - yi\~'"dxidyi < . 


0 


Xi - X2 
oPIa 


0 


Zl - 


yi - y2 

^Pla 


yi - X 2 


yi -y2 

j.PIOL 


Ala 


dyi 


\xi-yi\ ""dxidyi 
0(^l)dZl^ 0 / “ ^2 




j-Pjc 


dyi 


Now from the fact that c < a and c < b implies c < \/ab, the lemma follows. 

The following lemma is from [3, Theorem 3.5]. 

Lemma 7.4. Let Tn be the simplex defined in (5.5). Then for all h > —1, it holds that 


□ 


'T„(p 


[(^ ^n)An '5n—l) • • • ("52 '®l)] 


r(i + h)^ 


r(n(l + h) + 1) 


^n(l+h,) 


Proof of Theorem 3.14- Following the same notation and arguments as the proof of Theorem 
3.5 until (5.1), we have 


u! ||/n(') ll'H®" 

1 


< C- 




dsdr / dydzgn{s,yfi,x)gn{r,zfi,x)^K{iyi-Zi) 




i=l 


2 = 1 


Furthermore, by Cauchy-Schwarz inequality, we obtain 

« n 

dydz gn{s, y, f, x)gnir, z, t, x) JJ A{yi - Zi) 


/R2nd 


2=1 












< 


dydz gn{s, y, t, x)gn{s, z, t, x) JJ A{yi - Zi) 


i=l 

n 


1/2 


1/2 


x< / dydz gn{r,y,t,x)gn{r,z,t,x)\YKyi-^i) 


2=1 


Applying Lemma 7.3 to the above two integrals, we have 


where 

with 

Hence, 


'E2nd 


dydz gn{s, y, t, x)gn{r, z, t, x) JJ A{yi - Zi) < 


2=1 


n n 

2=1 2=1 


0 < s^(i) < Sa( 2 ) < ... < Sa{n) and 0 < rp(i) < rp( 2 ) < ... < rp(„). 


Qn !> n 

n\\\fn{-,-,t,x)\\‘^^n < / Jj7(si-’"i)( 0 (s) 0 (’^))^dsdr 

i=l 


< 


< 


n! 

J[0. 

^ a,l3,d,u,K 

1 

n\ 

2 j 


f 

n\ 

J[0, 

^ a,l3,d,i',K 

c? 


i=l 


]^7(si - Ti) {4>{sf^ + (j){rf^) dsdr 


i=i 


]^7(5i - ri)0(s)^®dsdr 


n\ 




(j){sy^ds 


_ 


I 0(s)^®ds 

T„it) 

r((20 + l)n + l) 


where Ct is defined in (2.2). Therefore, 


n!||/„(-,-,t,a;)||^®n < 


^ a,j3^d^iy,K.^t 

T{{2e + l)n + l) ’ 


and X]n>oconverges if 6* > —1/2. Finally, the condition 9 > —1/2, which 
is eqnivalent to k < 2q! — a/(3, guarantees both condition 6^ > —1 in Lemma 7.4 and the 
assumption k < 2a used in Lemma 7.3. This completes the proof of Theorem 3.14. □ 
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A Fox H-function 


Definition A.l. Let m,n,p, q be integers such that 0 < m < q,0 < n < p. Let a*, 6^ G C be 
complex numbers and let aj, (3j be positive numbers, i = 1,2,... ,p and j = 1,2,... ,q. Let 
the set of poles of the gamma functions r{bj + f3js) doesn’t intersect with that of the gamma 
functions r(l — a* — a^s), namely, 




_^-ai + k 

dik j k> 0,1, * • • 

CXi 


(A.l) 


for alH = 1, 2,..., p and j = 1, 2,..., g. Denote 

nr=i + ajs) nr=i r(l - - a^s) 




nr 


z=n+l 


T{aj + ais) ni=m+i r(l - - ajs) 


The Fox H-function 


TTm,n/...\ — TTm,n 
^p,q [^) = ^p,q 



(ai,Q;i) ■■ 

(dp, otp) 

z 

(&i,0i) •• 

■ {Kl^q) _ 


is defined by the following integral 


= ^ /_ n^:q^{s)z-^ds , zee. 


(A.2) 


(A.3) 


where an empty product in (A.3) means 1 and L in (A.3) is the infinite contour which 
separates all the points bji to the left and all the points Oik to the right of L. Moreover, L 
has one of the following forms: 

(1) L = L_ao is a left loop situated in a horizontal strip starting at point — cx) + i<pi and 
terminating at point —oo + z 02 for some —oo < 0i < 02 < exo; 

(2) L = L+oo is a right loop situated in a horizontal strip starting at point +cx) + and 
terminating at point oo + i 02 for some —oo < 0i < 02 < oo; 

(3) L = Li-foo is a contour starting at point 7 — ioo and terminating at point 7 + ioo for some 
7 G (—00, 00). 



Figure 1: Some illustrations of the path L. 

See Figure 1 for some illustrating paths. According to [22, Theorem 1.1], the integral 
(A.3) exists, for example, when 


q p 

A := > 0 and L = L_oo, (A.4) 

j=l i=l 
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or when 


n p m q 

a* := ^ ttj + ^ > 0 and L = Li^^. (A.5) 

2=1 2=n+l j=l j=m-\-l 

We call the poles 0 *^ (see (A.l)) of r(l — a* + ajs), z = 1,..., n, simple, if 

Vz,j, Q!j(l-ai + /c) 7 ^ Q!i(l-% + 05 j = 1,2,... ,n; /c, / = 0,1,.... (A.6) 

Similarly, the poles bji (see (A.l)) of T{bj + (3js), j = 1,... ,m, simple, if 

Wi,j, Pj{bi +k) ^ Pi{bj+ 1), i^j; i, j = 1,2,... ,m; k,l = 0,l,.... (A.7) 

Theorem A.2 (Part (hi) of Theorem 1.2 in [22]). Let L = Li^oo, a* > 0, z ^ 0, then 

is analytic on {z : jar^f^l < a*7r/2}. 


A.l Some elementary properties 

Property A.3 (Property 2.2 of [22]). If n > 1 and q > m, then 


h: 


m,n 




And similarly, if m > 1 and p > n, then 


h: 


m,n 


5*^1 ) 1 ,P—1 5 




Tjva-l.n I 


(Cti 


{bj,yj)2,q 


Property A.4 (Property 2.3 of [22]). It holds that 


TTm,n 

^p,q 


((2j ,0^4 ) 1 ,p 


_ TTn,m 


(1 bj,yj)l^q 

(1 


Property A.5 (Property 2.4 of [22]). For k > 0, it holds that 


H- 


m,n 

P^Q 


{ai,o>i)l^p 
{bj )l,g 


_ p Tjm,n 


(cij -jkot-i ) 1 ,p 
{bj ^kfSj)\^q 


A.2 Asymptotics at zero and infinity 


The following two theorems are some basic resnlts on the asymptotic expansions of the Fox 
H-fnnctions at zero and infinity. 

Theorem A.6. Suppose satisfies either A < 0, a* > 0 or A > 0. If the poles of 

P(1 — a, + Ojs) do not coincide, i.e., (A.6) holds, then when z —>■ oo, we have 


n oo 


i=l k=0 


(A.8) 


where 

hik 


■ 1 )^ nr.i r fi + |1 - g. + kffi) nUi*, r (i - - 11 - g. + <:]g) 
n?.„+i r {a, + |1 - a, + 4 ]g) n’™+i r (l - - |1 - a, + fc] a) 


k\ai 


(A.9) 
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Proof. See section 1.5 of [22] . (A.8) is asymptotic expansion 1.5.1 of [22]. □ 

Theorem A.7. Suppose satisfies either A < 0, a* > 0 or A > 0. When z ^ 0, we 

have the following two cases: 

(1) If the poles bji are simple (see (A.7)j, then 

m “ bj+i 

(A-io) 

i=i 1=0 


(2) If the poles bji (see (A.l ) ) ofT{bj + /3js), j = 1,... ,m, coincide: 

3i,j, I3fibi + k) = fifibj + 1), i ^ j; i, j = 1, 2,..., m; fc,/ = 0,l,---, 

then 


+ i\oez]\ 

’ ’ i=0 


N*,-l 

J* bi+l 


(A.ll) 


Here 


h* — 
IT-jl - 


(-1)' “ 1% + riLi^ 

nr.„+i r (oi - 16 ,+!||) n’™+. r (i - *,+ ib ,+!||) ’ 


(A.12) 


Yl'ji *'5 summation over j, I such that the bji do not coincide; Yl'ji summation over j, I 

such that bji coincide with order Nfi; 


n;,-i 


U* _ 

^jli - 




K, - 1 \ /n 


i I jl 


n 


m - 1 )' t 

where Nfi — 1 — n and n — i are orders of derivative; 


■ ]mhn}pr^-’''’mbii)r-'\ ( a . is ) 


Jn* 


3n* 


h;(s) = (s - bf J] r(6j + fts), m;(s) = (s - bf J] r(6j + (a.m) 


3=31 


3=31 


where b is the pole with order N*. 

Proof. See section 1.8 of [22] . (A.10) is asymptotic expansion 1.8.1 of [22]. (A.ll) is asymp¬ 
totic expansion 1.8.2 of [22]. □ 


A.3 A convolution property 

The following theorem is a simplihed version of Theorems 2.9 and 2.10 in [23], which is 
sufficient for our use. We need another parameter for the Fox H-function (A.2): 

P = ^bj-^ai + ^^-^. (A.15) 

j=l i=l 
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Theorem A.8. Let (a^, Ai,/ii) and (a^, A 2 ,/i 2 ) be the constants (a*, A,/i) defined in (A.5) 
(A.4) and (A. 15) for the following two Fox H-functions: 


Tjm.n 


X 




and 


H 


M,N 

P,Q 


X 


{di,5i)i^P \ 
(cj>7j)i,Q/ ’ 


respectively. Denote 


Ai = 
A 2 = 


mm 

l<i<n 


3fJ(aj) 


og 

5?(c,) 


mm 

l<j<M 7^- 


R - min 

, JJi — mm —-— 
l<i<m fij 

B 2 = min — 

i<i<iv Si 


iR{di 


with the convention that mm((;/)) = +cx). If either of the following four conditions holds 


(1) > 0 and > 0 and A 2 7^ 0; 


(2) > 0 and > 0 and Ai 7^ 0; 

(3) = Ai = 0, 3f?(/ii) < —1 and > 0; 

(4) 02 = A2 = 0 , 3f?(/X2) < —1 and al > 0 ; 

(5) = Al = 0, 3fJ(/ii) < —1 and = A2 = 0, ^(^ 12 ) < —I, 

and if 


Al + Bi > 0, A 2 + B 2 > 0, Al + A 2 > 0, Bi + B 2 > 0, 

then, for all z > 0, x E M, 


TTm+M,n+N 
^p+P,q+Q ' 


(q,^,Q: j)l^rii (^^i>Cki)n+l,p 

{bj ^/3j)i^rn, {bj ■,^j)m + l,q 


I 


{bj,gj)i,q 


H 


M,N ( f_ 


P,Q 


(diA)i,p\ dt 
(c:/>7i)i.Q/ t 


Proof. By Property A.4, 


H 


M,N (^ 


P,Q 


t 


(di,5i)q,p 

('=j,7j)l,Q 


ttNM 

^Q,P 


X 


(l-Cj, 7 j)l,Q 


(A.16) 


(A.17) 


If both Oi > 0 and > 0, then one can apply Theorem 2.9 of [23] with 7 = 0, a = 1, 
w = 1/x, and with the following replacements: N —>■ M, M ^ N, P ^ Q, Q ^ P, 
Cj —j- 1 — Cj, dj —)■ 1 — di. When al > 0 but 02 = 0 and A 2 7 ^ 0 , because both 2 ; and x are 
real numbers, Theorem 2.9 of [22] still holds (see the paragraph before Theorem 2.10 of [22] 
on p.59). Therefore, we have proved the theorem under conditions (1) or (2). 

If either of conditions (3)-(5) holds, we apply Theorem 2.10 of [23] in the same way. Note 
that the parameters fi for both Fox H-functions in (A.17) are equal. □ 
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A.4 Some integral transforms 

The following result on the Laplace transform of the Fox H-function is used in this paper. 
Denote 

POO 

(L/)(t):=/ fix)e-^^dx, t e C. 

Jo 

Theorem A.9 (Part of Corollary 2.3.1 of [22]). Suppose that a* > 0. Assume that w E C, 
a > 0 and a > 0 are such that 


a mm 


Then 


(ax^ 



{bj ,l3j)l^q) _ 


T<f<m fj. 


{t) = 


+ 3fJ(w) > —1. 


(A.18) 


-h: 








We also need the following convolution result, which is related to the Hankel transform 
of the Fox H-function. 

Theorem A. 10. (Corollary 2.5.1 of [22]) Let a* > 0 or a* = A = 0 and < —1. Let 
ri,w E C, r > 0 and a > 0 be such that 

+ ^(w) + T min > _i 

i<j<m /3j 

. 1 - 3?(ai) a 

T mm - > — — -I- 1, 

l<i<n Olj 2 


and 

Then for all a > 0 and b > 0, 

rixtrJ,iaixtr)H;^^^ (bf 


3?(r/) > -1/2. 


2ax V a 


w + 1 


T^m,n+l 

-”p+2,q 


X' 


2 
a 


dt 


rjo 


(bj,0j)l,q 


, X > 0. 


A.5 Riemann-Liouville fractional derivatives 


Recall that the Riemann-Liouville fractional derivative of order a E (0,1) is dehned in (4.1). 
The following results are used several times in this paper, which is a special case of part (i) 
of Theorem 2.8 of [22]. 

Theorem A.11. Let a E (0,1), ro G C and a > 0. If a* > 0 and (A.18) holds, then 


r^a Tjm.n / j.o 
-^0+^ ^ 


{bj ,0j)l^qj _ 


(x) = 


(bjAj)l,q,(-w+a,cr)J ■ 
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A.6 Some special cases 

The Mittag-Leffler function is a special case of the Fox H-function (see (2.9.27) of [22]): 




(A.19) 


The cos(-) function can be represented as (see [22, (2.9.8)]): 


cos(z) = ^/^ H 



(A.20) 


The Bessel function of the hrst kind Jn{z) is equal to (see (2.9.18) of [22]) 



(A.21) 


Another special case is (see (2.9.4) of [22]) 



(A.22) 


Acknowledgements 

Le Chen thanks Francesco Mainardi for some useful discussion and thanks Anatoly N. 

Kochubei for pointing out Pskhu’s paper [31]. 

References 

[1] R. Balan, D. Conus. Intermittency for the wave and heat equations with fractional noise 
in time. Ann. Probah. to appear (2015). 

[2] R. Balan, D. Conus. A note on intermittency for the fractional heat equation. Statist. 
Probab. Lett. 95 (2014), 6-14. 

[3] R. Balan, C. Tudor. Stochastic heat equation with multiplicative fractional-colored 
noise. J. Theor. Probab. 23 (2010), 834-870. 

[4] L. Chen. Nonlinear stochastic time-fractional diffusion equations on M: moments. Holder 
regularity and intermittency. Preprint arXiv:1410.1911, (2014). 

[5] L. Chen, R. C. Dalang. Moment bounds and asymptotics for the stochastic wave equa¬ 
tion. Stochastic Process. Appl. 125 (2015), no. 4, 1605-1628. 

[6] L. Chen, R. C. Dalang. Moments, intermittency, and growth indices for the nonlinear 
fractional stochastic heat equation. Stock. Partial Differ. Equ. Anal. Comput., to appear. 


2015. 


[7] L. Chen, R. C. Dalang. Moments and growth indices for nonlinear stochastic heat 
equation with rough initial conditions. Ann. Probab., (to appear), 2014. 


35 








[8] L. Chen, Y. Hu, D. Nualart. Nonlinear stochastic time-fractional slow and fast diffusion 
equations. Preprint, 2015. 

[9] Z. Q. Chen, K. H. Kim, P. Kim. Fractional time stochastic partial differential equations. 
Stochastic Process. Appl. 125 (2015), no. 4, 1470-1499. 

[10] R. C. Dalang. Extending the martingale measure stochastic integral with applications 
to spatially homogeneous s.p.d.e.’s. Electron. J. Probab., 4 (1999), no. 6. 

[11] K. Diethelm. The analysis of fractional differential equations. An application-oriented 
exposition using differential operators of Caputo type. Lecture Notes in Mathematics, 
2004. Springer-Verlag, Berlin, 2010. 

[12] S. Eidelman, A. Kochubei. Cauchy problem for fractional diffusion equations. J. Differ¬ 
ential Equations 199 (2004), no. 2, 211-255. 

[13] M. Foondun and D. Khoshnevisan. Intermittence and nonlinear parabolic stochastic 
partial differential equations. Electron. J. Probab. 14 (2009), 548-568. 

[14] M. Foondun, E. Nane. Asymptotic properties of some space-time fractional stochastic 
equations. Preprint at arXiv:1505.04615, (2015). 

[15] Y. Hu. Heat equations with fractional white noise potentials. Appl. Math. Optim. 43 
(2001), no. 3, 221-243. 

[16] Y. Hu. Chaos expansion of heat equations with white noise potentials. Potential Anal. 
16 (2002), no. 1, 45-66. 

[17] G. Hu and Y. Hu. Fractional diffusion in Gaussian noisy environment. Mathematics 3 
(2015), 131-152; doi;10.3390/math3020131. 

[18] Y. Hu, J. Huang, D. Nualart, S. Tindel. Stochastic heat equations with general multi¬ 
plicative Gaussian noises: Holder continuity and intermittency. Electron. J. Probab. 20, 
(2015), no. 55, 1-50. 

[19] Y. Hu, J. Huang, K. Le, D. Nualart, S. Tindel. Stochastic heat equation with rough 
dependence in space. Preprint arXiv:1505.04 924, 2015. 

[20] Y. Hu, J. Yan. Wick calculus for nonlinear Gaussian functionals. Acta Math. Appl. Sin. 
Engl. Ser. 25 (2009), 399-414. 

[21] D. Khoshnevisan, Y. Xiao. Harmonic analysis of additive Levy processes. Probab. Theory 
Related Pields 145 (2009), no. 3-4, 459-515. 

[22] A. A. Kilbas. and M. Saigo. H-transforms: theory and applications. Analytical Methods 
and Special Functions, 9. Ghapman & Hall/GRG, Boca Raton, FL, 2004. 

[23] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo. Theory and applications of fractional 
differential equations. North-Holland Mathematics Studies, 204. Elsevier Science B.V., 
Amsterdam, 2006. 


36 



[24] A. Kochubei. Diffusion of fractional order. Differential Equations 26 (1990), no. 4, 
485-492. 

[25] F. Mainardi, Y. Luchko, G. Pagnini. The fundamental solution of the space-time frac¬ 
tional diffusion equation. Fract. Calc. Appl. Anal. 4 (2001), no. 2, 153-192. 

[26] J. Mijena, E. Nane. Intermittence and time fractional stochastic partial differential 
equations. Preprint arXiv:1409.7468 (2014). 

[27] J. Mijena, E. Nane. Space-time fractional stochastic partial differential equations. 
Stochastic Process. Appl. to appear (2015). 

[28] D. Nualart. The Malliavin calculus and related topics. Second edition. Probability and 
its Applications (New York). Springer-Verlag, Berlin. 

[29] F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, editors. NIST handbook of 
mathematical functions. U.S. Department of Commerce National Institute of Standards 
and Technology, Washington, DC, 2010. 

[30] 1. Podlubny. Fractional differential equations, volume 198 of Mathematics in Science 
and Engineering. Academic Press Inc., San Diego, CA, 1999. 

[31] A. V. Pskhu. The fundamental solution of a diffusion-wave equation of fractional order. 
Izv. Math. 73 (2009), no. 2, 351-392. 

[32] W. Schneider. Completely monotone generalized Mittag-Leffler functions. Exposition. 
Math. 14 (1996), no. 1, 3-16. 

[33] W. Schneider, W. Wyss. Fractional diffusion and wave equations. J. Math. Phys. 30 
(1989), no. 1, 134-144. 

[34] J. Song. On a class of stochastic partial differential equations. Preprint 
arXiv:1503.06525v2, (2015). 

[35] E. Stein. Singular integrals and differentiability properties of functions. Princeton Math¬ 
ematical Series, No. 30 Princeton University Press, Princeton, N. J., 1970. 

[36] E. Stein, G. Weiss. Introduction to Fourier analysis on Euclidean spaces. Princeton 
Mathematical Series, No. 32. Princeton University Press, Princeton, N. J., 1971. 

[37] V. V. Uchaikin, V. M. Zolotarev. Chance and stability. Stable distributions and their 
applications. Modern Probability and Statistics. VSP, Utrecht, 1999. 

[38] D. V. Widder. The Laplace transform. Princeton Mathematical Series, v. 6. Princeton 
University Press, Princeton, N. J., 1941. 


37 



